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ABSTRACT 

Context. An important benchmark for current observational techniques and theoretical modeling of exoplanet’s atmosphere is the 
transit of Venus (ToV). Stellar activity and, in particular, convection-related surface structures, potentially cause fluctuations that can 
affect the transit light curves. Surface convection simulations can help the interpretation of ToV and also other transits outside our 
Solar System. 

Aims. We used realistic three-dimensional (3D) radiative hydrodynamical (RHD) simulation of the Sun from the SiAGGER-grid and 
synthetic images computed with the radiative transfer code Optim3D to provide predictions for the transit of Venus (ToV) in 2004 
observed by the satellite ACRIMSAT. 

Methods. We computed intensity maps from RHD simulation of the Sun and produced synthetic stellar disk image as a observer 
would see, accounting for the centre-to-limb variations. The contribution of the solar granulation has been considered during the ToV. 
We computed the light curve and compared it to the ACRIMSAT observations and also to the light curves obtained with solar surface 
representations carried out using radial profiles with different limb-darkening laws. We also applied the same spherical tile imaging 
method used for RHD simulation to the observations of center-to-limb Sun granulation with HINODE. 

Results. We managed to explain ACRIMSAT observations of 2004 ToV and showed that the granulation pattern causes fluctuations 
in the transit light curve. We compared different limb-darkening models to the RHD simulation and evaluated the contribution of 
the granulation to the ToV. We showed that the granulation pattern can partially explain the observed discrepancies between models 
and data. Moreover, we found that the overall agreement between real and RHD solar granulation is good, either in term of depth or 
Ingress/Egress slopes of the transit curve. This confirms that the limb-darkening and the granulation pattern simulated in 3D RHD 
Sun represent well what is imaged by HINODE. In the end, we found that the Venus’s aureole contribution during ToV is ~ 10~ 6 
times less intense than the solar photosphere, and thus, accurate measurements of this phenomena are extremely challenging. 
Conclusions. The prospects for planet detection and characterization with transiting methods are excellent with access to large a 
amount of data for stars. Being able to explain consistently the data of 2004 ToV is a new step forward for 3D RHD simulations that 
are becoming essential for the detection and characterization of exoplanets. They show that the granulation have to be considered as 
an intrinsic incertitude, due to the stellar variability, on precise measurements of exoplanet transits of, most likely, planets with small 
diameters. In this context, the importance of a comprehensive knowledge of the host star, including the detailed study of the stellar 
surface convection is determinant. 

Key words. Planet-star interactions - Sun: granulation - Techniques: photometric - stars: atmospheres - hydrodynamics - radiative 
transfer 


1. Introduction 


The detection of exoplanet orbiting around stars is now a routine 
with more that 1500 confirmed planets and about 3400 candi - 
dates (as in fall 2014 from http://exoplanets.org, I Wright etaD 
201 il l. The study of exoplanet atmospheres can generally be 
classified into two main categories: (i) the measurement of time- 
varying signals, which include transits, secondary eclipses, and 
phase curves, and (ii) the spatially resolved imaging and spec¬ 
troscopy. These techniques allow to explore the property of ex¬ 
oplanet’s atmospheres and understand the diversity of chemi¬ 
cal compositions of exoplanets, their atmospheric processes, in¬ 


ternal structures, and formation conditions (Madhusu dhan et alJ 
l20T4h . 

A transit event occurs when an (or more) exoplanet move 
across the face of host star, hiding a small part of it, as seen 
by an observer at some particular vantage point. The exoplanet 
blocks some of the starlight during the transit and creates a peri¬ 
odic dip in the brightness of the star. The amount of light reduc¬ 
tion is typically ~ 1%, 0.1%, and, 0.01% for Jupiter-, Neptune - 
and Earth-like planets, respectively (iBorucki & Summers! 19841) . 
This depends on the star-planet size ratio and on the duration 
of the transit, which, in turn, depends on the planet’s distance 
from the host star and on the stellar mass. This method is mostly 
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sensitive to large, short-period exoplanets, which display the 
largest atmospheric signatures and make them popul ar targets 
for ti me-varying characterization studies (Ma dhusudhan et all 
120141) . During both the primary and secondary transits, the 
planet-star radius ratio can be deduced, and, for sufficiently 
bright stars, also the exoplanet’s mass can be determined 
from the host star’s radial velocity semi-amplitude dMislis et al.l 
120121) . If the mass and/or radius of the exoplanet is good enough, 
also the atmospheric density can be predicted giving hints to the 
current understanding of exoplanet formation processes. Even 
more information can be retrieved by measuring the depth of the 
transit: the radiation coming from the host star passes through 
the exoplanet’s atmosphere and it is absorbed to different de¬ 
grees at different wavelengths. This wavelength-dependent depth 
of the transmission spectrum is directly linked to the absorp¬ 
tion features imprinte d on starlight transmitted through the ex- 
oplanet’s atmosphere dSeager & Sasselovl l2000: Brown 200 lh . 
However, the transit method has several drawbacks: (i) to be ob¬ 
served, the exoplanet must cross directly the line-of-sight from 
the Earth (i.e., the star-planet system must be edge-on with re¬ 
spect to the observer) and this is the case only for a small mi¬ 
nority of exoplanets, (ii) Another problem concerns the inter¬ 
transit duration that can last months to years, reducing drasti¬ 
cally the number of detectable exoplanets, (iii) There is also 
a chance to observe the alignment with a background eclips¬ 
ing binary that cause blends difficult to disentangle dTorres et al] 
2011). (iv) A last potential complication is the host star sur¬ 
face inhomogeneities caused by the presence of stellar granula¬ 
tion (e.gJChiavassa et al.ll2Q14h or the magnetic starspots (e.g., 
IStrassmeie3l20Q9l). 

A very important benchmark for current observational tech¬ 
niques and theoretical modeling of exoplanet’s atmosphere is the 
transit of Venus (ToV). Venus is in many ways an Earth analogue 
and it’s the closest planet transit we can see in our Solar System. 
ToV events are rare, as they occur in pairs 8 years apart, each pair 
separated by 121.5 or 105.5 years, alternating between descend¬ 
ing node (June pairs: 1761/1769, 2004/2012,...), and ascending 
node (December pairs: 1631/1639, 1874/1882, 2117/ 2125,...). 
The 2004 and 2012 transits of Venus represent thus a crucial 
leap on human-life timescale, as the modern imaging technolo¬ 
gies available allow for the first time a quantitative analysis of 
the atmosp heric phenomena of the planet associated to its transit 
dTanga et al.ll2012l) . Spatial- and ground-ba sed data have been 
collected and studied by different authors. iTanga et all (120121) 
performed the first photometric analysis of the "aureole'Q us¬ 
ing various ground-based system s and extracted f undam ental pa¬ 
rameters of Venus’s atmosphere. [Garcia Munoz & Millsl ( 2012 ) 
compared images of Venus during Ingress/Egress to 2004 data 
and provided guidelines for the investigation of the planet’s 
upper haze from vertically-u nresolved photometric measure¬ 
ments. lEhrenreich et al.l ( 2012 ) provided a theoretical transmis¬ 
sion spectrum of the atmosphere of Venus dominated by car¬ 
bon dioxide absorption and droplets of sulfuric acid composing 
an upper haze layer (features not expected for a Earth-like ex¬ 
oplanet) that could be tested with spectroscopic observations of 
2012 transit. 

In this work, we present the transit light curve predictions 
obtained from 3D surface convection simulation of the Sun and 
compare them to the ToV occurred in 2004. Furthermore, we 
compare 3D predictions to limb-darkening models of the Sun, 


1 At the Ingress and Egress of ToV, the portion of the planet’s disk 
outside the Solar photosphere has been repeatedly perceived as outlined 
by a thin and bright arc, called "aureole" 


largely used in the planetary community, as well as to the real 
observations of the granulation pattern on the Sun. The aim of 
this work is to test, whether the 3D surface convection approach 
to the transit technique gives reliable results for the benchmark 
case of the ToV, which also is important for the interpretation of 
other transits outside our Solar System. 


2. From the hydrodynamical simulation of the Sun 
to spherical tile imaging 

The Sun is a natural reference for studying other stars and so¬ 
lar models are essential for many studies in stellar and plan¬ 
etary astrophysics. Many of the observable phenomena occur¬ 
ring on the surface of the Sun are intimately linked to convec¬ 
tion. Moreover, the atmospheric temperature stratification in the 
optical thin region, where the emerging flux form, is also af¬ 
fected by the interaction between radiative and convective en¬ 
ergy transport. To account for all these aspects, it is impor¬ 
tant to use realistic 3D radiative hydrodynamical (RHP) sim¬ 
ulations of convection jst ein & Nordlund 1998; iFrevtag et al.l 
l2012t ICarlsson et al.l 120041: IVogler et al.l 120041) . RHD simula¬ 
tions of the Sun are widely used and compared in detail with 
observations by the astrophysical commun i ty using different nu- 
merical cod es (e.g.. iNordlund et al.ll20Q9b lAsplund et al.ll2009l 
ICaffau et al .11201 iVlP ereira et al. II 201 ih. 


2. 1 . Stellar model atmospheres 


We used the solar simulation from the STAGGER-grid of real¬ 
istic 3D RHD sim ulations of stellar convection for cool stars 
(iMagic et al.1120131) . This grid is computed using STAGGER-code 
(originally developed by Nordlund & Galsgaard 19950, and con¬ 
tinuously improved over the years by its user community). In 
a cartesian box located at the optical surface (i.e., r ~ 1), the 
code solves the time-dependent equations for conservation of 
mass, momentum, and energy coupled to an realistic treatment of 
the radiative transfer. The simulation domains are chosen large 
enough to cover at least ten pressure scale heights vertically and 
to allow for about ten granules to develop at the surface, more¬ 
over there are periodic boundary conditions horizontally and 
open boundaries vertically. At the bottom of the simulation, the 
inflows have a constant entropy and pressure. The simulations 
employ realistic input physics: (i) the equation of state is an up¬ 
dated version of the one bv lMihalas et al. ( 19881) : (ii) the radia¬ 
tive transfer is calculated for a larg e number over wavelength 
point s merged into 12 opacity bins (lNordhmdl [l982: Ska rtlienl 
l2000l) with opacities in cluding cont i nuous absorption and scat- 

1(l2010l) and line opacities as 

20Q8), 

1200 ll) of atomic lines. The 


tering coefficients from Hayek et a 


described from Gustafsson et al. 
VALD-2 database dStempels et al 


abundances employed in the computation of the Sun simulation 
are the latest chemical composition bv lAsplund et al.l ( 2009 ). 


2.2. Three-dimensional radiative transfer 

The computation of the monochromatic emerging intensity is 
performed using the 3D pure local thermal eq uilibrium radiative 
transfer code Optim3D dChiavassa et al.ll2009l) and the snapshots 
of the RHD simulation of the Sun (see Table [T]). The radiative 
transfer equation is solved monochromatically using pretabu¬ 
lated extinction coefficients as a function of temperature, density, 

2 http://www.astro.ku.dk/~kg/Papers/MHD_code.ps.gz 


Article number, page 2 ofiTTI 













































































A. Chiavassa et al.: A new view on exoplanet transits: Transit of Venus described using three-dimensional solar atmosphere 
Table 1. 3D simulation of the Sun from STAGGER-grid used in this work. 


< T eff >^_ 

[Fe/H] 

log,? 

x, y, z-dimensions 

x, y, z-resolution 

M* 

R* 

Number of tiles ^ 

[K] 


[cgs] 

[Mm] 

[grid points] 

[M g ] 

[R©] 

over the diameter 

5768.51 (Sun) 

0.0 

4.4 

7.7x7.7x5.2 

240x240x240 

1.0 

1.0 

286 


a Horizontally and temporal average of the emergent effective temperatures from [Magic et al J ( 20131 ) 

b Mile = — ;- Re ■ 

U1C x,y -dimension 


and wavelength based on the same extensive atomic and molec¬ 
ular continuum a nd line opacity data as t he latest generation of 
MARCS models dGustafsson et al.ll2008h . We assume zero mi¬ 
croturbulence and model the non-thermal Doppler broadening 
of spectral lines using only the self-consistent velocity fields 
present in the 3D simulations. We employed th e sam e chemical 
composition of the 3D RHD simulations (lAsplund et al .11 11200 9). 
The wavelength interval computed is 6684.0 ± 0J_ A. The de¬ 
tailed methods used in the code are explained in lChiavassa et al.l 
(120091) . 

2.3. Spherical tile imaging 

The computational domain of each simulation represents only 
a small portion of the stellar surface (see, e.g., the central 
panel of Fig. [TO] and Fig 1 in IChiavassa et al.l 20121) . To ob¬ 
tain an image of the whole solar disk accounting for limb- 
darkened effects , we em ploy the same tili ng method expla i ned in 
Chiavassa et al.l (120101) and used also in IChiavassa et al] (120121 
20141) . Optim3D has been used to compute intensity maps from 
the solar simulation for different inclinations with respect to 
the vertical, p= cos(0)=[l.OOOO, 0.9890,0.9780, 0.9460, 0.9130, 
0.8610, 0.8090,0.7390,0.6690, 0.5840,0.5000,0.4040,0.3090, 
0.2060, 0.1040], with 6 being the angle with respect to the line- 
of-sight (vertical axis), and for 30 representative snapshots of 
the simulation adequately spaced apart so as to capture several 
convective turnovers. 

We then used the synthetic images, chosen randomly among 
the snapshots in the time-series, to map onto spherical surfaces 
accounting for distortions especially at high latitudes and lon¬ 
gitudes cropping the square-shaped intensity maps when defin¬ 
ing the spherical tiles. The computed value of the 6 -angle used 
to generate each map depended on the position (longitude and 
latitude) of the tile on the sphere and was linearly interpolated 
among the inclination angles. The random choice of snapshots 
avoid the assumption of periodic boundary conditions and the 
resulting tiled spherical surface displays an artifactual periodic 
granulation pattern (Fig. [TJ. The total number of tiles (Vtiie) 
needed to cover half a circumference from side to side on the 
sphere is Mile = ^ y _dtoLion = 286 ’ where R ° ( transformed in 
Mm) and the x, y -dimension come from Table [TJ We generated 
50 different synthetic solar disk images. Since the number of 
tiles necessary to cover the sphere is larger than the total number 
of representative snapshots of RHD simulation, some tiles may 
randomly appears several times on the solar disk at different in¬ 
clinations. Therefore, the 50 synthetic images are not completely 
independent, but we assume that they are a good enough sta¬ 
tistical representation to estimate the granulation changes with 
respect to time during the ToV. Fig. [5] displays the fluctuation 
of the granulation structures for a particular cut in the synthetic 
disk images. The solar disk intensity fluctuates at the scale of 
the tile and the maximum intensity varies between [3.74-3.87] 
xlO 6 erg cm -2 s -1 A -1 . 



X [Rsun] 

Fig. 1 . Synthetic solar disk image computed at 6684.0 + 0.1 A of 
the RHD simulation of Table Q] The intensity range is [0.0-3.79 x 
10 6 erg cm -2 s -1 A -1 ]. We generated 50 different synthetic solar disk im¬ 
ages to account for granulation changes with respect to time. 

3. Modeling the transit of Venus of 2004 

June 8th 2004 was a very important and rare date for monitor¬ 
ing the first of two close by (200 4 and 2012) transits of Venus 
visible from Earth, before 2117. [Schneider et all ( 2006 ). here¬ 
after SPW2006, took the chance and obtained space-based solar 
irradiance measurements with ACRIMSAT spacecraft. 

3.1. Observations with ACRIMSAT 

The satellite mission ACRIMSAT^, launched in 1999 and 
founded by NASA through the Earth Science Programs Office 
at Goddard Space Flight Center, measured 0.0961% reduction 
in the solar intensity caused by the shadow of the Venus during 
his transit in 2004 (SPW2006). These measurements were taken 
with the Active Cavity Radiometer Irradiance Monitor (ACRIM 
3) instrument ( IWillson & Mordvinov1200 3). ACRIM 3 was de¬ 
signed to provide accurate, highly precise, and traceable radiom- 
etry over decadal timescales, to detect changes in the total energy 
received from the Sun by the Earth. ACRIM 3 provides 32 second 
exposure shutter measurements characterized by open (observa¬ 
tions) and closed (calibration) measurements. The measurement 
precision is 0.01% (SPW2006). 


3 http://acrim.jpl.nasa.gov 
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Fig. 2. Cut at x=0, y>0 for 50 different synthetic solar disk images (one 
example in reported Fig. [I]). The intensity profiles are normalized to the 
intensity at the disk center (R=0.0). The green line is the RHD <3D> 
average profile. 


ACRIMSAT cannot observe the Sun continuously because 
the line-of-sight is occulted by the Earth during its orbit around 
the Sun. In SPW2006, these periodic interruptions were deter¬ 
mined for time intervals spanning the ToV using a definitive 
a posteriori orbital ephemeris derived from a contemporaneous 
epochal satellite element set provided by the North American 
Aerospace Defense Command. SPW2006 reports (Figure 2 in 
SPW2006) the apparent path of the center of Venus as seen from 
ACRIMSAT with respect to the heliocenter as the planet crossed 
the face of the Sun. ACRIMSAT s observations are interrupted 
for ~30 minutes by Earth occultation during each of its ~100 
minute orbits. The duration of 2004 ToV is ~5.5 hours. 

In this work, we use the trajectory coordinates of Venus extracted 
from SPW2006 to simulate the passage of the planet on the syn¬ 
thetic solar disk of Fig. [T] 



3.2. Evaluation of Venus’s aureole with HINODE 

Telescope images, with instrument typically up to 15 to 20 cm, 
report cusp extension tending to transform the thin crescent of 
Venus into a ring of light since the mid 1 8th ce ntury (for detailed 
historical reviews see lLinkll 1 9691: lEdsonll 1 963h . While these phe¬ 
nomena can be ascribed to light diffusion, during transits the 
contribution of refraction (3-4 orders of magnitude stronger) 
shows up. Traditionally, its discovery is attrib uted to Lomonosov 
(who lived between 1711 and 1765; Marov 2005 ). As the ToV 
are rare, also the data concerning the aureole are sparse. How¬ 
ever, thanks to the coordinated giant leap in the observation of 
2004 ToV, as well as the modern technologies allowed a quan¬ 
titative analy sis o f the a tmospheric phenomena associated to the 
ToV (iTanga et al.1120121) . The authors performed for the first pho¬ 
tometric analysis and obtained spatially resolved data providing 
measurements of the aureole flux as a function of Venus’s lati¬ 
tude along the limb. 

The aureole is a distorted (strongly flattened) image of a 
portion of the solar photosphere. As such the surface brightness 
of the photosphere is the same in the aureole. For this reason, 
although the aureole is very thin (< 40 km above the cloud 
deck) it is clear the overall contribution to the transit photometry 
deserves a specific evaluation. To do this, we needed high spatial 
resolution images of Venus during a transit, and, only HINODE 
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Fig. 3. Top panel: Image (1024 x 1024 pixels) of Venus (dark circu¬ 
lar subaperture) transiting the Sun (brighter areas) in June 5th 2012 as 
seen by HINODE SOT/SP spectropolarimeter (total intensity Stoke pa¬ 
rameter /) in the red filter centered at 6684.0 + 0.1 A (Red filter). The 
intensity range is 0.-2.7 x 10 3 DNpixel -1 s -1 (where DN is data num¬ 
bers, iLites et aLll2013h . The intensity plotted is the square root of the 
intensity to better display the "aureole". Central panel: Same figure as 
above with the contrast highly increased to fit the radius of the planet, 
yellow circle, with the procedure described in the text (R fit = 274.91 
pixels). The center of the Sun is located in South-East direction. Bot¬ 
tom panel: Three dimensional view of the images above. The aureole 
contribution is visible in the region just outside the solar disk. 
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spacecraft could provide that. Unfortunately, the satellite was 
launched after 2004 and thus we used observations carried out 
during ToV of 2012, assuming that there were no changes in 
Venus’s atmospherebetween2004 and2012. 

HINODE dKosugi et al.l 120071) is a joint mission between the 
space agencies of Japan, United States, Europe, and United 
Kingdom launche d in 2006, and since then, the spe ctropo- 
larimeter SOT/SP dLites et al.ll200lb iTsuneta et al.ll2008l) allows 
the solar community to perform spectropolarimetry of the solar 
photosphere under extremely stable conditions. SOT/SP pro¬ 
vides 0.32 arcsec angular resolution measurements of the vector 
magnetic fields on the Sun through precise spectropolarimetry 
of solar spectral lines with a spatial resolution of 0.32 arcsec 
and photometric accuracy of 10 -3 (ilchimoto et al.llll2008l) . 

We obtained HINODE observations of June 5th 2012 ToV 
from through the HINODE data on-line archive. We used 
SOT/SP Stokes I images (i.e., the total intensity measured) 
in the red filter (6684.0 ±0.1 A). The data are reduced and 
calibrated for all the slit positions (i.e., corresponding to "levell" 
data). The solar diameter in that period was 31.7'(~951.1"for 
the radius), the geometrical distance at the Sun surface, seen 
under one arcsecond was /?©soho/951.1 x 0.109 = 732.1 km 
(with /?©soh o = 696342 km measured by SOHO spacecraft, 
lEmilio etafl 120121) . The spatial resolution of the image is 
1024x1024 pixels (Fig. [3]) and 1 pixel = 0.109", the size of each 
image is then 111.61 x 111.61" or 81.71 x 81.71 Mm, with 1 
pixel equal to 79.80 km. 

The contribution from the aureole was computed as follows: 

1. The dark circular subaperture corresponding to Venus in 
Fig. [3] ( centr al panel) was fitted using Newton-based fit 
method ( Prat t 1987 ). We selected the points where approx- 
imatively resides Venus’s radius (yellow numbered crosses 
in the figure). The fitting method returned the center of the 
planet as well as its radius (or impact parameter), a in the 
figure, with a Venus radius 7?fi t =274.91 pixels. 

2. The position of the center and the radius known, we extracted 
a ring of 40 pixels in the outer regions of Venus within the 
region ± 20 pixels. 

3. We extracted 215 radial intensity profiles within the ring re¬ 
gion (a representative example is reported in Fig. [4]). 

4. We fitted the aureole contribution using a gaussian law (3 
free parameters) for the peak and a straight line equation 
(2 free parameters) for the background. Then we integrated 
over the radius the intensity subtended by the gaussian fit¬ 
ting curve for each of the 215 radial profiles. We obtained 
the radial contribution of the aureole. 

5. We repeated this procedure for 50 other images during the 
Ingress and Egress of ToV with similar results. 

Fig. © shows the typical signal of the aureole in the ring 
region selected close to Venus’s radius. The peak is signifi¬ 
cant and much larger than the background signal (always lower 
than ~20 DN pixel -1 s -1 ). The contribution of the aureole cor¬ 
responding to the 215 radial intensity profiles is in the range 
^aureole = [0-4 - 6.0] x 10 5 DN pixel -1 s -1 with an average value 
of < /aureole >= 1.2 x 10 5 DNpixel -1 s -1 (Table©. 

To relate the aureole intensity to the full Sun, we use the 3D syn¬ 
thetic solar disk image of Fig. © First, we define /mnode as the 
total emerging intensity from Fig [3] (top panel), / syn thgiobai as the 
total emerging intensity from the synthetic Sun of Fig. © and 
4ynthcut as intensity emerging from a cut in the synthetic disk of 



Fig. 4. Representative radial intensity profiles within the ring region 
(R fit +20m pixels) in the image of Fig.© which includes the aureole con¬ 
tribution. The blue dotted curve corresponds to the observations, while 
the green solid line to the best fit (see text for details). The vertical red 
dashed line is the location of Venus’s radius. 

Fig. ©corresponding to the same size and shape of HINODE ob¬ 
servation of Fig © Then the relative intensity due to the aureole 

IS /ratio = /aureole X (/synthcut//synthglobal) X (1/ /HinodeX knowing that 

(/synthcut//synthgiobai) = 0.0118. Table ©reports values of the order 
of 10 -6 , far too weak to be detected, also considering the preci¬ 
sion level of ACRIMSAT (~ 10 -5 ) and the height of the transit 
(10 -4 , Fig. [6]). We therefore do not consider the aureole in our 
modeling. 

Table 2. Integrated intensity contributions (/ aU reoie) of maximum, min¬ 
imum and average values from the aureole (Fig. 0 , together with the 
related ratio (/ ra tio) to the total emerging intensity from the synthetic 
Sun of Fig.© 



^aureole X 1 (P 

[DN pixel” 1 s” 1 ] 

/ratioXlO” 6 

Average 

1.2 

1.0 

Max 

6.0 

5.4 

Min 

0.4 

0.3 


3.3. Modeling the transit with 3D RHD simulation of the Sun 

The observations give the transit parameters such as: the geomet¬ 
rically derived transit depth and contact times, the ratio between 
Venus’s and Sun’s radii, and the positions of Venus with respect 
to the heliocenter as a functions of time given the spacecraft and 
planetary orbits. In the following, we explain how we modeled 
the ToV using RHD simulation of the Sun as a background. The 
steps are: 

1. we use the synthetic solar image (Fig.© as the background 
emitting source; 

2. Venus is modeled by a disk of radius 7?v enus = 0.0307/?© ~ 
8.5 pixels (Fig.© with total intensity 0.0961% of / syn thgiobai 
(percentage flux decrement during 2004 transit, SPW2006); 

3. the planet disk crosses the solar disk following the apparent 
trajectory of Venus as seen from ACRIMSAT (see SPW2006 
for more details and Fig.©: 

4. the intensity of the system is collected for every transiting 
step. 
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Fig. 5. Top panel: Synthetic solar disk image in the visible of the RHD 
simulation with the different positions of ToV as seen from ACRIMSAT. 
The computed wavelength is the same as in Fig.[T] The unusual apparent 
trajectory of Venus is induced by the spacecraft orbit (see text). Bottom 
panel: enlargement from top panel. 

The intensity map of Fig. [5] is a statistical representation of 
the Sun at a particular moment during the transit. The apparent 
trajectory of Venus on the Sun is characterized by (see SPW2006 
for more details): (i) a periodic spatial and temporal modulations 
in the location of Venus as it traverses the solar disk caused by 
the projection on the Sun of the shifting of the line-of-sight in¬ 
duced by the spacecraft orbit; (ii) a vertical (north-south) ampli¬ 
tude variations resulting from the near-polar orbit of the space¬ 
craft; (iii) an horizontal (east-west) component manifesting in 
nonlinear spacings in the planetary position along its projected 
path in equal time intervals due to orbit plane not being on the 
line-of-sight direction to the Sun. 

The temporal fluctuations is taken into account using the 50 
different synthetic solar disk images computed. The resulting 
light curve from the ToV is reported in Fig. [6j The light curve 
has been averaged over 50 possible realizations of solar disk. The 
decrease of the Sun’s light is due to the planet occultation and de¬ 
pends on the ratio between the respective size of Venus and the 
Sun as well as on the brightness contrast. The magnitude of this 
intensity variation depends on the ratio of the size of Venus on 
the brightness contrast between the solar disk and the planet. As 




Universal time (Hours from 08 June 2004 Oh UT) 

Fig. 6. Top panel: ToV 2004 light curve (black dashed line) for the 3D 
RHD simulation of the Sun (Fig. [5]) compared to the photometric ob¬ 
servations (red dots with error bars, [Schneider et al.l120060 taken with 
ACRIM 3 mounted on ACRIMSAT. The observati ons have been normal¬ 
ized to 1, dividing the original data (Figure 3 of [Schneider et al.1120061) 
by 1365.88 W/m 2 . The light curve has been averaged (with the green 
shade denoting maximum and minimum values) over 50 different syn¬ 
thetic solar disk images to account for granulation changes with respect 
to time. Bottom panel: same as above but for one particular realization 
among the 50 synthetic solar disk images. 


a result of the reflective spacecraft parallactic motion of Venus, 
the planet’s apparent path crosses regions of the Sun with sur¬ 
face brightnesses very close to the stellar limb. This affects the 
slope of the light curve that is more shallow during Egress with 
respect to the Ingress. This aspect is very constraining for transit 
modeling and crucial to validate the RHD simulation used. 
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Fig. 7. Cu t at x=0, y>0 in different synthetic solar disk images either 
built from IClaretl (120001) (light blue curve), IMagic et aD (120141) (red 
curve), RHD <3D> profile (black curve, from Fi g |2]) and HINODE 
data for the quiet sun (green curve, used in Section [331 . The intensity 
profiles are normalized to the intensity at the disk center (R=0). 

The data obtained in 2004 for the ToV (see Section 13.11) 
has been provided to us by Glenn H. Schneider (University of 
Arizona, private communication). Unfortunately, 2012 ToV data 
are not helpful in our work because they had too large error 
bars. Venus occulted 0.0961% of the total area of the Sun pho¬ 
tosphere, and, ACRIM 3 provided calibrated observations with 
gaps resulting from the Earth occultations during the satellite 
orbit. Fig. [6] shows that the agreement between our RHD simu¬ 
lation and the data is very good. The root-mean-square deviation 

(RMS D = ^ f=1 ^ 1 ^ ; where y* are the predicted values, yt the 
observed values, and n the number of measurements) varies be¬ 
tween [9.1-9.5]xl0 -5 for all the 50 possible realizations of solar 
granulation and is equal to 9.1xl0 -5 for the average profile plot¬ 
ted in Fig. [6] At the bottom of the light curve (~8.4h UT), there 
are some aperiodic variations. SPW2006 reported that this may 
result from intrinsic changes in the total solar irradiance over the 
same time interval or may arise as Venus occults isolated regions 
of the photosphere differing in local surface brightness (e.g., 
sunspots or smaller spatial scale localized features). Indeed, it 
is po ssible that magnetic starspots c ontaminate the transit signal 
(e.g.. ISilvall2QQ3b iFanza et al.ll2009l) . In the end, accounting for 
the granulation fluctuations allows to reduce the scatter on the 
observed points at about 7h, 9h, 1 lh UT (Fig. [6l bottom panel). 
Our approach is statistical and aims not to fully explain the ob¬ 
served data but to show that our RHD simulation of the Sun is 
adapted (in terms of limb-darkening and emerging flux) to in¬ 
terpret such a data and that the granulation fluctuations have an 
important effect of the light curves during the transit. The evalu¬ 
ation of the granulation contribution is done in Section [3~4l 

3.4. 3D RHD granulation versus limb-darkening models 

In the literature, in general, the contribution of the star to the 
transit is approximated by a parametric representation of the ra¬ 
dial limb-darkened surface profile. This approach neglects the 
surface inhomogenieties caused by the granulation pattern, the 
intensity brightness decrease is taken into account, depending on 
the limb-darkening law used. As a consequence, the slope of the 
light curve in the Ingress and Egress may be affected differently. 
In SPW2006, the authors used this method and examined the 


causes that could affect their deficiency in light curve fitting of 
the transit data. In their work, they argued that instrumental mea¬ 
surement errors and intrinsic solar variations compared to imper¬ 
fections in the limb-darkening model itself may be the cause. In 
this context, RHD simulation of the Sun brings a new insight on 
the impact of the granulation thanks to the excellent matching 
showed in Fig. [6] The numerical box of a RHD simulation is ro¬ 
tated with a determined mu-angle (see Section [23]) before com¬ 
puting the emerging intensity. There is not a pre-adopted limb 
darkening function and the final result is the direct output of the 
emerging intensity at a particular mu-angle. 

Table 3. Limb-darkening coefficients used in this work with the law 
Ifilh = 1 - z£= 1 °k (l - /i k/2 ^J ( Claretll2000i) and a star with solar param¬ 
eters or for the <3D> profile of the RHD simulation (Fig. [2] green line) 
or the intensity profile of HINODE data (Fig. [7] green line). 


Reference 

a\ 

Cl2 

a 3 

a/i 

Claret et al. 2000 

0.5169 

-0.0211 

0.6944 

-0.3892 

Magic et al. 2014 

0.5850 

-0.1023 

0.5371 

-0.2580 

Fit RHD <3D> 

-3.8699 

13.3547 

-15.4237 

6.1910 

Fit HINODE data 

18.32 

-38.49 

36.05 

-12.20 


We tested different limb-darkening laws represent the Sun. 
For this purpose, we built synthetic images representing the so¬ 
lar surface using_radial profiles obtained with the limb-darkening 
law of [Claret 2000 (based on ID model atmospheres of lKuruca 
19791 ): I^/h = 1 - i a k (l _ y^ /2 )> expressed as the variation 
in intensity with //-angle that is normalized to the disk-center 
(Ijj/h)- The limb-darkening coefficients for a star with the stel¬ 
lar parameters of the Sun and for the visible range are in Table [3] 
The reader should note that [Magic et afl (12014l) used the same 
RHD simulation of Table [T] to extract their coefficients. More¬ 
over, we used the law of lClaretl to fit the average RHD profile 
< 3D > of Fig. [2] (green curve) and reported the coefficients in 
Table [3] 

The normalized intensity profiles of the different laws are re¬ 
ported in Fig. [7] 

For the circular subaperture representing the disk of Venus 
we used a constant flux. Fig. [8] ( top) shows the three different 
limb-darkening m odels: iMagicTet all (120141) (RMSD = 9.0 x 
10 -5 , red curve), IClare if <l2000h (RMSD = 9.3 x lO" 5 , light 
blue curve), and the Claret law for the RHD <3 D > profile 
(RMSD = 9.3 x 10 -5 , purple curve). Fig. [8] (bottom) displays 
that the central depre ssion o f Venus’s transit obtained with the 
limb-darkened disk of Claret ( 2000 ) (light blue curve) is system¬ 
atically undersized with respect to the transit with the 3D RHD 
simulation (black curve), while there is a better agreement with 
IMagic et aD (12014l) model (red curve) and with the RHD <3 D > 
limb-darkening. Since the magnitude of the intensity variation in 
the light curve depends on the brightness contrast between the 
solar disk and Venus, it is noticeable that RHD synthetic disk 
image differs from limb-darkening models with particular em¬ 
phasis at the stellar limb. The contribution of the granulation to 
the ToV can be retrieved by the absolute difference between the 
limb-darkening of the RHD <3 D > profile and the full RHD 
ToV (purple curve in bottom panel of Fig. [8]). We computed the 
RMSm = 1.7 x 10 -5 of this difference and found that it is ~5 
smaller than the RMSDs of the different models. Even if the 
granulation contribution is not large enough to explain the whole 
observed discrepancies, it is source of an intrinsic noise due to 
the stellar variability, that may affect precise measurements of 
exoplanet transits of, most likely, planets with small diameters. 
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Universal time (Hours from 08 June 2004 Oh UT) 

Fig. 8. Top panel: Comparison of different ToV light curves. The black 
curve is the light curve combining the RHD simulation of the Sun 
(Fig. \5\ and Fig. (6]> . The red a nd li ght blu e curve s represent the limb- 
darkening Ma gicet ahl (120141) and lClaretl (120001) . respectively. In the 
end, the purple curve represents the limb-darkening carried out using 
Claret’s law for the RHD <3 D > profile of Fig. [2] (green curve). In all 
the cases, a constant value is used for Venus intensity. Bottom panel: 
Difference between the above curves with respect to the transit curve of 
the RHD simulation of Fig. [6] 

3.5. 3D RHD granulation versus observed solar granulation 

The aim of this Section is to create a solar disk image using the 
same procedure of RHD simulation (Section 12.31) but with, as 
input, the granulation pattern observed on the Sun. For this pur¬ 
pose, we recovered HINODE spectropolarimeter SOT/SP data 
product collected between December 19 and December 20 2007 
from High Altitude Observatory (HAO) websit^J The data are 
reduced and calibrated for all the slit positions (i.e., correspond¬ 
ing to "levell" data). We used SOT/SP Stokes I images (i.e., 
the total intensity measured) integrated at the continuum wave¬ 
lengths very close the Fe I line pair at 6300 A but avoiding the 
spectral lines: the emerging intensity is then compatible, in tex 
of wavelengths probed, to the intensity map of Fig. Q] and the 

4 http://www.csac.hao.ucar.edU/csac/archive.jsp#hinode 
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Fig. 9. Sun total intensity image (370 x 370 pixels) as seen by HINODE 
SOT/SP spectropolarimeter at different Ap and at continuum wave¬ 
lengths close to 6300 A. Top panel is the center of the Sun and bot¬ 
tom panel is the southern limb. The intensity range is, from the top, 
18-27 x 10 3 - 17-27 x 10 3 - 14-21 x 10 3 - 18-25 x 10 3 DN pixel" 1 s" 1 . 
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HINODE Red filter used in Fig. [3] The original images of 1024 
pixel-spectrograph slit parallel to the North-South polar axis of 
the Sun was successively located at 20 latitudes and scans of 
370 steps of 0.16" in the direction perpendicular to the slit were 
performed. At each position in the slit, the images have been in¬ 
tegrated for about 5 seconds and for a total time of about 30 min¬ 
utes (ie., 370 steps). This provides us with 20 (1028x370 pixels) 
images allowing a complete recording of the center-to-limb vari¬ 
ations of the solar disk. However, as the fields of view of the 
images are partly superimposed we extracted from the HINODE 
images a set of 17 (370x370 pixels) images forming a contin¬ 
uous South-North stripe in the solar southern hemisphere. The 
solar diameter in that period was 32.5'(~975.3"for the radius), 
the geometrical distance at the Sun surface, seen under one arc- 
second was Rosoho/975.3 = 713.9 km. Each pixel is 0.16", the 
size of each image (Fig. [9]) is then 59.20x59.20" or 42.18x42.18 
Mm, with 1 pixel equal to 1 14 km. 

However, the images were too large compared to the syn¬ 
thetic maps (Fig. [TO] central panel) we used for tiling the spher¬ 
ical surface of the Sun (Fig. [TJ with RHD simulation. We then 
divided each image by a factor 10 in the East-West axis and by a 
factor 10 in the North-South polar axis. We ultimately obtained 
several 37 x 37 images of ~4.2 Mm (Fig. [TCI top panel) in 83 p 
sub-steps. 

In the end, we applied the tiling procedure of Section [23] and ob¬ 
tained the ToV for a granulation pattern based on HINODE ob¬ 
servations (Fig. [TO] bottom panel). The deviation from ACRIM- 
SAT observations is RMSD = 8.7 x 10 -5 . Fig. ITTI displays the 
comparison between the light curve obtained from the RHD sim¬ 
ulation (Fig. [6] black dashed line) and the one calculated from 
the ToV of Fig.flOl (bottom panel). The absolute difference be¬ 
tween the ToV curves (Fig. eh bottom panel) shows a typi¬ 
cal value of ~ 5.0 x 10 -5 , which is close to what is reported 
in Fig. [8] (bottom). However, the agreement is better, in partic¬ 
ular for the Ingress/Egress slopes of the transit. This confirms 
that the limb-darkening and the granulation pattern simulated in 
the 3D RHD solar simulation represent well what is imaged by 
HINODE. This comparison with HINODE data is important be¬ 
cause allows another inspection of RHD simulation using data 
with smaller noise with respect to ToV data. 

Following what we did in Section [3~4l for the RHD <3D> pro¬ 
file, we used the limb-darkening law of Claret to fit the inten¬ 
sity profile obtained from HINODE data (Fig. [7] green curve). 
The coefficients are reported in Table [3j It is also noticeable in 
Fig.[7]the good agreement between the Hinode and the RHD pro¬ 
files. To obtain the contribution from the granulation, we com¬ 
puted the limb-darkening ToV and subtracted it to the ToV of 
the HINODE image. The resulting signature of the granulation 
has an RMS Hinode = 3.0 x 1CT 5 , which is comparable to the 
RMS in = 1.7 x 10 -5 due to the granulation signal in RHD ToV. 
Both values are however lower than or of about the same or¬ 
der of the smallest error bar of ACRIMSAT (2.7 x 10 -5 ), and 
thus the granulation signal cannot be unveiled with these data. 
Nevertheless, It should be note that there is a possibility that the 
small scale variability produced by the RHD simulation may be 
slightly underestimated due to the snapshot chosen. 




0 12 3 4 

x [Mm] 
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4. Conclusions and implications for transit studies 
of extrasolar planets 

Modeling the transit light curve of Venus as seen from ACRIM¬ 
SAT implies the importance to have a good representation of 
the background solar disk and of the planet itself. For the Sun, 


Fig. 10. Top panel: Cut of 37 x 37 pixels from Fig. [9] top panel. 
Central panel: intensity map of a small portion of the Sun surface 
at p=l from the RHD of Tabled] The intensity range is [2.0-4.3 x 
10 6 erg cm -2 s -1 A -1 ]. Bottom panel: ToV on the solar disk image con¬ 
structed with the tiling procedure of Section [23] but with real granula¬ 
tion observations (top panel and Fig. [9} of HINODE (see text for de¬ 
tails). A constant value is used for Venus intensity. 
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Fig. 11. Top panel: ToV light curves for RHD simulation (Fig. [5] and 
Fig. [6]) and HINODE SOT/SP spectropolarimeter data (Fig. [Toj bottom 
panel). A constant value is used for Venus intensity. Bottom panel : Dif¬ 
ference between the HINODE ToV with respect to the RHD simulation. 


we used the realistic, state-of-the-art, time-dependent, radiative- 
hydrodynamic stellar atmosphere of the Sun from the Stagger- 
grid. For Venus, we extracted and fitted radially averaged profiles 
from high spatial resolution images of HINODE. 

We provided synthetic solar disk image using t he sp herical 
tile imaging method already applied in IChiavassa et al.l (120141 
120121 12010) . We applied a statistical approach to show that 
our RHD simulation of the Sun is adapted (in terms of limb- 
darkening and emerging flux) to interpret such a data and that 
the granulation fluctuations have an important effect of the light 
curves during the transit. Our procedure successfully explain 
ACRIMSAT observations of 2004 ToV and showed that the gran¬ 
ulation pattern causes fluctuations in the transit light curve. We 
compared different limb-darkening models to our RHD simula¬ 
tion and evaluated the contribution of the granulation to the ToV. 
We showed that the granulation pattern can partially explain the 
observed discrepancies between models and data. However, the 
amplitude of the signal is close to the typical differences between 


the different models used and it is difficult to disentangle the sig¬ 
nal of the granulation with the data used in this work. Moreover, 
we found that the Venus’s aureole contribution during ToV is 
~ 10 -6 times less intense than the solar photosphere, and thus, 
accurate measurements of this phenomena are extremely chal¬ 
lenging. 

In the end, we applied the same spherical tile imaging proce¬ 
dure to the observations of center-to-limb Sun granulation with 
HINODE and built a solar disk image. We then compared the 
resulting transit curves and found that the overall agreement be¬ 
tween real and RHD solar granulation is very good, either in term 
of depth or Ingress/Egress slopes of the transit. This confirms 
that the limb-darkening and the granulation pattern simulated in 
3D RHD Sun represent well what is imaged by HINODE. 

3D RHD simulations are now established as realistic descrip¬ 
tions for the convective photospheres of various classes of stars. 
They have been recently employed: (i) in the prediction of differ¬ 
ential spectroscopy during exoplanet transits to reconstruct spec¬ 
tra of small stellar surface portions that successively become hid¬ 
den behind the planet (iDravins et al.ll2Ql4h : (ii) or to assess the 
transit light curves effect of the limb-darkening issued from 3D 
RHD simulations with respect to ID models. 

The ToV is an important benchmark for current theoretical mod¬ 
eling of exoplanet transits. Being able to explain consistently the 
data of 2004 transit is then a new step forward for 3D RHD sim¬ 
ulations. This implies that our tiling procedure can be adapted 
to exoplanet transits. Moreover, our statistical approach in this 
work allowed to show that the granulation have to be considered 
as an intrinsic incertitude, due to the stellar variability, on precise 
measurements of exoplanet transits of, most likely, planets with 
small diameters. A possible solution to reduce this incertitude 
is the repetition of the transit measurements several times rather 
than take a single snapshot. 


The prospects for planet detection and characterization with 
transiting methods are excellent with access to a large amount 
of data for different kind of stars either with ground-based 
telescope s suc h as HATNet(Hungarian Automated Telescope 
Networ k. iBakos et al.l 120041) . NGTS (Next Generation Transit 
Survey. Iwheatlev et~al.ll2013l) . TRAP PIST (TRAnsiting Planets 
and Planeteslmals Small Telesc ope. iJehin et al. 1201 lh . WASP 
(Wide Angle Search for Planets, iPollacco et alJl2QQ6b or space- 
base d missions like CHEOPS (CHaracte rizing ExOPlanet Satel¬ 
lite, iBroeg et aPl2013l:lkoch et al.ll2010h. COROT (C onvection, 


Rotation and planet ary Transits. [Baglin et al.ll2006bllah . Kepler 
dBorucki et al.l2010h . TESS (Transit Exoplanet Survey Satellite, 
iRicker et all201Ch. PLATO (PLAnetary Transits and Oscillation 
of stars, iRauer et al.l2014l) . 

In this context, 3D RHD simulations are essential for a detailed 
quantitative analysis of the transits. Indeed, the interpretation 
and the study ofthe imp act o f stellar granulation is not limited to 
the Sun (iNordlund et al.l2009h but RHD simulations cover a sub¬ 
stantial portion of the Hertzsprung-Russell diagram (Magic et all 
120131: iTrampedach et al.l T2013I: iLudwig etak 12009 ). including 
the evolutionary phases from the main sequence over the turnoff 
up to the red-giant branch for low-mass stars. 
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